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文章 编号 :1005-3085(2010)04-0669-10 


二 阶 椭圆 问题 的 最 小 二 乘 扩 展 混合 有 限 元 方法 * 


T 胜 ， 陈 焕 贞 1 
(山东 师范 大 学 数学 科学 学 院 ， 济 南 250014) 


H E 为 克服 最 小 二 乘 混合 元 方法 在 数值 模拟 具 小 扩散 系数 或 低 渗透 率 问题 时 ， 应 对 扩散 系数 求 逆 带 来 
的 困难 ， 本 文 基于 最 小 二 乘 与 扩展 混合 元 的 思想 ， 对 一 类 刻画 扩 获 、 渗 透 过 程 的 二 阶 椭圆 问题 建 
立 了 最 小 二 乘 扩展 混合 元 格式 ， 证 明了 格式 的 稳定 性 和 收敛 性 质 。 论 证 表明 ， 该 格式 具有 无 需 对 
小 扩散 系数 求 送 ， 较 好 的 克服 了 小 扩散 系数 带 来 的 困难 ， 能 同时 高 精度 通 近 未 知 函 数 ， 梯 度 及 其 
通 量 ， 有 限 元 空间 无 需 满足 LBB 条 件 ， 刚 度 和 矩阵 对 称 正 定 等 最 小 二 乘 方法 和 扩展 混合 元 方法 的 
良好 性 质 。 数 值 算 例 说 明了 所 提 算 法 的 有 效 性 。 
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HN CR (mn = 1,2,3) 为 有 界 凸 域 ,边界 T= TpUTw (Tp # 9) Lipschitz ER, ARE FI! 
二 阶 椭圆 边 值 问题 
—V -(AVu(2)) = f(z), TER, 
u(x) = 0, zerp, (1) 
(—AVu(z))-n =0, ZETN， 


其 中 Vu 表示 函数 以 的 梯度 ，V .gq 表示 向 量 函 数 g 的 散 度 ，n 是 的 单位 外 法 向 量 ， 甜 阵 A = 
(aij(x))P 21, TE 人 对 称 正定 ， 且 存在 正常 数 a, 和 a*， 使 得 


aX X S XTAX Lax X% YxER”. (2) 


对 上 述 二 阶 椭圆 问题 人 们 提出 了 很 多 数值 模拟 方法 ， 如 有 限 元 方法 由 ， 混 合 有 限 元 方 
法 回 ， 扩 展 混 合 有 限 元 方法 B4， 最 小 二 乘 有 限 元 方法 和 最 小 二 乘 混合 有 限 元 方法 6? 等 。 扩 展 
混合 有 限 元 方法 由 于 具有 能 同时 高 精度 逼近 未 知 函 数 ， 未 知 函 数 的 梯度 和 流体 的 通 量 ， 适 合 处 
理 具有 复杂 边界 条 件 和 小 系数 问题 的 优点 ， 因 此 在 地 下 水 动力 学 问题 的 模拟 中 得 到 了 广泛 的 应 
用 。 但 是 该 方法 仍然 要 求 有 限 元 空间 满足 LBB 相 容 性 条 件 ， 限 制 了 空间 的 选择 ， 并 且 所 形成 
的 方程 组 是 非 正 定 的 ， 计 算 复杂 。 为 更 好 地 数值 模拟 实际 的 地 下 水 动力 问题 ， 人 们 希望 构造 一 
种 既 能 具备 扩展 混合 有 限 元 方法 的 优点 又 不 需要 满足 LBB 相 容 性 条 件 ， 计 算 简 单 的 高 效 数 值 
模拟 方法 。 
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本 文 将 最 小 二 乘 思想 和 扩展 混合 有 限 元 方法 相 结合 ， 提 出 并 分 析 了 最 小 二 乘 扩展 混合 有 
限 元 方法 。 该 方法 继承 了 最 小 二 乘 和 扩展 混合 有 限 元 方法 的 优点 ， 即 格式 能 同时 高 精度 逼近 
未 知 函数 ， 未 知 函数 的 梯度 以 及 流体 的 通 量 ， 所 形成 的 有 限 元 方程 组 是 对 称 正定 的 且 条 件数 
为 O(h-?)， 易 于 计算 ， 有 限 元 空间 无 须 满足 LBB 相 容 性 条 件 ， 构 造 更 灵活 。 数 值 算 例 说 明了 


该 方法 的 可 行 性 和 有 效 性 。 


2 ERER 
首先 ， 给 出 几 个 用 到 的 空间 的 定义 
A= (LQ), V= {ve HQ) :v=0 Tp E}, 
H(div;Q) = {q € (L7(Q))": V -q € L2(0))}, 
W = {q € H(div;Q):q-n=0 Æ Tn E}. 
| H llog = (4, 人 二 H E€ A, 
lo lhe= (lv la+ IV la)? vev, 
la laanm= (1a a+ V- ala)? qew, 
其 中 (.,-) 表示 空间 LO) R (L) PHAR. 
记 入 = 一 Vu, o = Ad, WW (1) 可 变形 为 关于 w AM o 的 一 阶 方程 组 
AA—o=0, LER, 
A+Vu=0, «EQ, 
V-o-f=0, TEQ, 
u=0, 2 E 工 D， 


o-n=0, rETN. 


一 阶 形式 (3) 的 最 小 二 乘 泛 函 可 定义 为 


J(v, ma) =(V-q—-F,V-9-f) + (ut Vu, u+ Vv) + (Ap — q, Ap — 9). 


与 (3) 等 价 的 最 小 二 乘 极 小 问题 为 : 求 (w 和,o) EV x A x W， 使 得 


J(u, à o) = oe a J(v, p,q). 
相应 的 变 分 形式 为 : K (u, 和,o) EV x A x W， 使 得 
a(u, A, 030, 4,q) = (f,V- 49), Vue, Be A, qe W, 


其 中 


a(u,r,0;0, 1,9) = (V -0, V - q) + (A+ Vu, u + Vv) + (AA — o, Ap — q). 


(5) 


(6) 
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为 证 明 变 分 形式 (5) 的 解 是 存在 唯一 的 ， 我 们 首先 说 明 双 线性 形式 a(; .) 的 强制 性 。 
定理 1 双 线 性 形式 -) 是 强制 的 ， 即 存在 常数 C > 0， 使 得 对 任意 的 ve V, KEA, qE 
W， 有 
a(v, h, q; V, H, q) > c( | v l?a + I H lêo + il q lraiv:n) ). (7) 


证 明 由 (6) 式 ， 知 
a(v, 1,950,H,9) = (V -4,V -q) + (u+ Vv, u + Vu) + (Au — q, Au — q) 
= (V -4, V - q) + 28(q, Vv) + 67(v, v) — 26(Vv, q) — B7(v,v) 
+(H, u) + 2(u, Vv) + (Vu, Vo) + (Ap, Au) — 2(Au,a) + (3,9), (8) 
其 中 6 > 0 是 适当 选取 的 常数 。 应 用 Green 公式 ， 并 进一步 整理 ， 有 
a(v, 1,950, p,q) = (V-q— pv, V -q — Bu) + (4,9) — 2(Ap + BV2,q) 
+(Ap + BVv, Au + BVv) — (Ap + BV», Au + BVv) 
+(Ap, Ap) + 2(Vv, p) + (u, u) + (Vv, Vu) ~ 6?(v,v) 
= (V-q— bv, V -q — Bv) + (Au + Bu ~q, Au + Bv ~ q) 
+(u + (E — BA)Vv, u + (E — BA)Vv) + 26(AVv, Vv) 


—P2(AVv, AVv) — B? (Vv, Vu) — 6? (u,v) 


> 26(AV», Vu) — B?(AVv, AVv) — 6?(Vu, Vv) — B?(v, v), (9) 
SOP E ERMER. 
由 于 v(x) = 0, z ETp， 所 以 由 Poincaré-Friedrichs 不 等 式 知 ， 存 在 正常 数 Cr， 使 得 
lv lõa < Cr Il Vo lõ, - (10) 
又 由 (2) 知 
2B(AVv, Vv) > 28a.(Vv, Vv), (11) 
B?(AVu, AVv) < B?a*? (Vv, Vv). (12) 
从 而 由 (10),(11) #1 (12), 
a(v, u, q; v, 1,9) > B(2a, — B(1 + Cr + a**)) |] Vo HBa- (13) 


$ B=a,/(1+Cpr +a**), 则 有 


a,? 


一 *2 = e—a 
B(2a, — B(1+ Cr +a*°)) 1G; pa" > 0， (14) 


因此 


aU; p, gv, pq) Z C || Vu lia > C Ilo lita - (15) 
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Xha) 的 定义 (6)， 显 然 有 


a(v, p,q;v, 4,9) > |V -a llo (16) 
a(v, p,q; v, p,q) > | p+ Vv lao, (17) 
a(v, Hq; V, H, q) 2 | Ap -q lba s (18) 


从 而 由 (15) 和 (17)， 有 


ellie = I e+ Vv- Vv én 


IA 


2 I H + Vv én +2 | Vu ao < Calv, p, q; VU, H, q), (19) 
进一步 由 (18) 和 (19)， 有 


lellaao = I g- Aut An len 


IA 


2 || Au 一 ao +2 I| Ax lon < Calv, p, g; v, p,q). (20) 


综合 (15), (16), (19) 和 (20) 便 证 得 (7) 成 立 。 
定理 2 WE 72(Q)， 则 问题 (5) 存在 唯一 解 (w Na) EV X Ax W. 
证 明 考虑 空间 (VA,W)， 其 范 数 定义 为 


fv let te loe + lla lava, VEV, nehA, geWw. 


由 定理 1 知 双 线 性 形式 a(.;.) 在 空间 (V, A, W) 中 是 强制 的 ， 而 且 由 Cauchy PERDE a(-;-) 在 
(V,A,W) 中 是 有 界 的 。 男 外 ， 当 f € LCR, REER = (f,:) 也 是 有 界 的 ， 因 此 
由 Lax-Milgram 引 理 知 (5) 的 解 是 存在 唯一 的 。 

定理 3 设 六 72(Q)， 则 格式 (5) 是 稳定 的 。 

证 明 在 (5) 中 令 v=w, w=A, q=0, Wheel, E 


(Ilfsu IIo + 1A Neg + Ilo Waive) 


Calu, r,0;u,A,0) = (f, V-o) 


IA 


< C| F lõa tel V-o laos C I F l,o (21) 


A 


从 而 格式 (5) 的 稳定 性 得 证 。 
3 ”最 小 二 乘 扩展 混合 有 限 元 逼近 
为 简单 起 见 ， 我 们 假定 Q 是 R? A SCRE 中 的 凸 多 面体 。 有 是 Q 的 一 个 剖 分 ， 


即 
X= |] K, 


KETh 
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Eh = max{diam(K) : K € Ty}. Va, An MW, DREV, AMW 的 有 限 维 子 空间 ， 具 有 下 
面 的 逼近 性 质 


inf {| v — vn lon +h I| v — vr la} < Chow [le+1,0, (22) 
VhEVh 

inf || u- La llog < CRY || p |Ir+1,0, (23) 
HhEAnh 

inf {|i g- an llog +h || @— an llagan } < CR Ila lst, (24) 
qnEWh 


其 中 ,r,s 为 非 负 整 数 ，v e H*+1(Q), we (At (O))", q € (H'O). Va, A 和 Wi 的 一 
个 标准 的 选 法 是 k, + 和 s 次 的 分 片 多 项 式 空 间 ， 即 


Vi = {un € C°(Q) : un [KE Pe(K), VK € Th, vh=0 在 Tp E}, (25) 
An = {pn € (C°(Q))” : un |KE (Pr(K))", YK € Th}, (26) 
Wn = {an € H(div;Q) : gn |x € (P3(K))", VK € Th, Hen =O ETN E}, (27) 


其 中 Pn,(K), m= k,7,s 是 单元 KK 上 的 m 次 多 项 式 空间 。 
定义 (5) 式 的 有 限 元 逼近 格式 为 : 求 (un, An, on) E Vn x An x Wh， 使 得 


aluh, Ap oh; Vh, Ma, qh) = (f,V cah) Vurn€ Vn, wane An, qn € Wh. (28) 


AV, CV, An CA, Wh CW， 所 以 由 定理 1 知 问题 (28) 的 解 存在 唯一 。 
下 面 ， 我 们 将 检验 由 格式 (28) 生成 的 线性 方程 组 的 条 件数 。 为 此 ， 假 定 剖 分 二 是 正则 的 且 
满足 道 假定 ， 即 存在 常数 5 > 0， 使 得 


h 
二 一 一 KET. 
diam(K) 一 Oy: See ae (29) 


在 假定 (29) 下 ， 一 般 的 标准 有 限 元 空间 都 满足 下 面 的 逆 性 质 
lunli < Ch? || vn llon, Yun E Vh, (30) 
IV -arloa < Ch || gn llon, Van € Wa. (31) 


Bobir 0} A {by} BE ETAL Va, An MW BOSE, WUE RE AY (un, Anon) € 
Va x Ah Xx Wi, A 


L M N 
th = Somer, M=} Gb =) GY. 
l=1 i=1 j=l 


记 In|, 上 和 | 分 别 表示 向 量 (m, set NL); (1, The , EM) 和 (G1, Ck , CN) 的 L? 范 数 。 在 假定 (29) 
下 ， 常 见 的 有 限 元 空间 都 满足 如 下 的 假定 : PEERK o, bi 和 (i = 1,2)， 使 得 


oih"|n| < |] vp llo < a2h”|n|, Vur € Va, (32) 
Bik” El < |l pn lloa < Boh |E|, V Hn E An, (33) 


wh" IC < ll an lo < yeh" |Cl, Van E Wh. (34) 
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定理 4 假定 (30)-(34) 成 立 ， 则 线性 方程 组 (28) 的 条 件数 为 O(P2)。 
证 明 由 定理 1，o(.;) 的 有 界 性 以 及 (30) 和 (31)， 知 
Ci( Ilon loa +I a lõa + Ia lêo ) 
< a(Yh, Hh, gp; Vh, Hh, gh) < Coh-*( |] va lia + |] Ha lotlar lba). 
再 结合 (32), (33) 和 (34), A 
Crh” (In? + IEP? + II?) < alun, Has ghi Vh, pn gn) < C2h®™? (lm 十 长 2 二 |， 


此 即 证 得 (28) 的 条 件数 为 O(h2)。 
进一步 ， 我 们 得 出 下 面 的 误差 估计 。 
定理 5 tue Ht*+t1(Q), Ae (HHN), o € (Hs+1(Q))*， 则 有 


|lu — unr llao + |] A— An lon + lo ~ on ycaivo) 
< C(h* || u ller AHE |] A lri +° || o lls+rnn )- 
证 明 由 (5) 和 (28)， 可 得 到 误差 方程 


a(u — Un, À — Àh, — Ohi Uh, ay Gh) = 0, Von © Vi, th © An, qh © Wh. 
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(35) 


(36) 


Bur E Vn, Ar CA, Moy CW, DH yu, 入 和 cz 的 标准 有 限 元 插值 ， 则 有 下 面 的 插值 估 


itl 
|| u — ur llog +h || u — ur lao < Ch*** || u Ilevt,0, 
| A— Az llo < CREE 站 入 ao， 
lc 一 cr lloa +A || o ~ or laivo) < ChY lo lley - 
由 定理 1 和 (36)， 有 


lur = ur Mig + Il Ar Arlao+lon 一 careaivo) 





< Ca(up — Ur, Àn — ÀL, Oh — Or Un — Ur, Àh — ÀI Oh — O1) 
= Ca(u — ur, À — Mo — oF; Un — UT, Àh — ÀL, on — OF) 
< C(u- ur ĝa + IA- Ar Ba +o- or ravim)? 

x (|| ua — ur llia + Il An — àz lba + |l on 一 cz ean 

进一步 ， 利 用 (37), (38) 和 (39)， 可 得 
| ur — Ur lia + |] An — Az llo + I| on — or lyeaiv;o) 
< C( || u- ur lio + |] A— Ar llog + || o- oF laaivo)) 
< C (hë | u letro HAF |] A lleia +h? |] o lleria). 


由 (40) 并 利用 三 角 不 等 式 便 可 证 得 (35)。 


(37) 
(38) 


(39) 


(40) 
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4 ”数值 算 例 
为 验证 所 提 格 式 的 有 效 性 ， 我 们 首先 在 区 域 [0, 1] 上 考虑 一 维 问题 
| 一 (auz)z = cm TE (0, 1), (41) 


u(0)=0, u(1)=0, 


其 中 
= a f= mee + 32). 
容易 验证 方程 的 真 解 为 
peered. Aa Pee 6 2 eed), 


100 


取 An 为 分 片 常数 多 项 式 空间 ，Vi AW), 为 分 片 线性 多 项 式 空间 ， 用 Matlab 编程 计算 近似 
fi Un, An 和 on， 得 出 真 解 和 近似 解 的 误差 以 及 最 小 二 乘 扩展 混合 有 限 元 格式 的 收敛 阶 ， 详 见 



























































表 1 和 表 2。 
Rl: 误差 
HEA 
0.01 0.005 0.0025 0.00125 
{| u — un ljg 0.024040625 | 0.012039385 0.006024433 | 0.003013399 
| \— An lloo 0.024040621 0.012039385 0.006024433 | 0.003013399 
{| ec — on mcaivio) 0.000620657 | 0.000310923 0.000155610 | 0.000077842 





















|| u — un llig 
| A — An llo,n 


lc — on adivoa) 















通过 表 中 数据 可 以 看 出 ， 当 步 长 h 折 半 递 减 时 ， 误 差 也 折 半 递减 ， 误 差 的 收敛 阶 达到 了 最 
优 ， 这 说 明 本 文 所 提 格 式 是 有 效 的 。 进 一 步 ， 我 们 绘 出 真 解 和 最 小 二 乘 扩展 混合 有 限 元 解 ( 取 


步 长 h = 0.01) 的 比较 图 ， 见 图 1 至 图 3。 
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图 1: u5 un 的 比较 图 图 2: 入 与 和 的 比较 图 
0.015 = 
h=0.01 
0.01 pane J 
R, 
0.005 
ot wS 
N 
-0.005 F Ni 1 
N 
-0.01 X, 
N 
-0.015} \ 
\ 
-0.02} SN 4 
sa = \ 
0.0257} 一 A J 
-0.03 , ， 
0 0.2 04 0.6 08 1 








图 3: oS op, 的 比较 图 
EKO = [0, 1] x [0, 1] 上 考虑 下 列 二 维 问题 


| -V-(aVu)=f, TER, (42) 
u(x) = 0, ZE On, 


其 中 
1 1 1 1\2 1\? 
u= 16(1 — 2)ey(1— 9) (5 Ta arctan ( 20( = — (z— 5) = (v- 5) D) a=3. 
H Matlab FEW EEM un, An 和 cm， 得 出 真 解 和 近似 解 的 误差 以 及 最 小 二 乘 扩 展 混合 
有 限 元 格式 的 收敛 阶 ， 详 见 表 3。 
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表 3: || u—ua loa, IA- Arlon llo 一 ohllo,n 的 误差 及 收敛 阶 








节点 数 








通过 表 中 数据 可 以 看 出 ， 误 差 的 收敛 阶 达到 了 最 优 ， 与 理论 结果 一 致 。 进 一 步 ， 我 们 绘 出 
最 小 二 乘 扩 展 混 合 有 限 元 解 un, A 与 真 解 入 的 比较 图 ， 见 图 4 至 图 5。 





图 4: 与 wi 的 比较 图 ( 左 图 为 真 解 ， 右 图 为 数值 解 ) 





图 5: A oH An, on 的 比较 图 (上 面 两 图 分 别 为 真 解 和 , og， 下 面 两 图 分 别 为 数值 解 An, on) 
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Least-squares Expanded Mixed Finite Element Methods for 
Second-order Elliptic Problems 


DING Sheng, CHEN Huan-zhent 


(School of Mathematical Sciences, Shandong Normal University, Jinan 250014) 


Abstract: To overcome the difficulties when calculating the inverse of a small diffusive coefficient 
when simulating the diffusive problems within a low permeability zone by a least-squares mixed finite 
element method, we use an expanded mixed finite element instead of a mixed finite element and de- 
velop a least-squares expanded mixed finite method for second-order elliptic problems. We prove the 
stability and convergence for the proposed procedure. The analysis shows that the method inherits 
the advantages of least-squares and expanded mixed finite element methods, such as it is not necessary 
to calculate the inverse of a small diffusive coefficient; approximating the unknown, its gradient and 
its flux directly; the finite element space being free of LBB condition as well as the stiff matrix being 
symmetry and positive-definite. Numerical tests are performed to confirm the theoretical analysis. 
Keywords: second-order elliptic problem; least-squares method; expanded mixed finite element 
method; optimal-order convergence 
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